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Abstract— In this paper, ft technique is described which uses multiple 
surface acoustic wave (SAW) devices in parallel to reduce the acquisition 
time of a direct sequence spread-spectrum communication system. 
Analysis of system performance in both the search and lock modes is 
presented, and key quantities such as probability of false alarm, 
probability of correct detection, mean dwell time, and mean time to lose 
lock are derived. 

I. Introduction 

IN this paper, a rapid acquisition technique, based upon the 
use of parallel processing the received waveform of a direct 
sequence spread-spectrum communication system with multi- 
ple surface acoustic wave (SAW) convolvers, is presented. 
Because the processing is indeed done in parallel, a decrease 
in acquisition time from a comparable serial search scheme is 
achieved. 

A derivation of the probabilities of false alarm and correct 
detection for both the search mode and the lock mode are 
presented, as well as a Markov chain analysis of an overall 
search-lock strategy. Expressions are derived for such key 
system .parameters as probability of entering lock and average 
time to loss of lock. Finally, some numerical results are pre- 
sented, showing the variation of false alarm and detection 
probabilities v/ith the number of parallel convolvers and with 
the ratio of received energy per bit to noise spectral density. 

The paper is divided into six sections, with a description of 
the system presented in the next section, the derivations of 
probability of false alarm and probability of correct detection 
presented in Section III, a Markov chain analysis of the search 
and lock modes* behavior given in Section IV, and numerical 
results given in Section V. Conclusions are then presented 
in the last section. 

II. System Description 
The acquisition scheme can be best described by referring 
to Fig. 1. There are two modes of operation, a search mode 
and a lock mode. Consider the search mode first, during which 
all switches are in position 1. (Ignore, for the moment, the 
recirculating delay line circuit at the input to the system.) 
Notice that there are N convolvers, each of duration 2T 
seconds, where it is assumed that M chips of the spreading 
sequence span T seconds. The full period of the spreading 
sequence, L t is divided into subsequences each of length Af, 
and for simplicity, it is assumed that L/M is an integer. 
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Each .of the N convolvers has as a reference input one of the 
subsequences of length M. Initially, let us assume that the 
total phase uncertainty at the start of acquisition, say K 
chips, is spanned by the MN phase positions of the N con- 
volvers. (Notice that even if L is much greater than MN, 
this uncertainty condition could still be satisfied.) This 
guarantees that the correct phase position of the received 
waveform is "seen" by one of the N convolvers somewhere 
in the integration interval of the convolvers. 

The convolver outputs are sampled at times t = T + fT c , 
where T = MT C and / = 0, 1, M - 1. The length of each 
convolver is chosen to be 2T. With these parameter values, 
if the N reference inputs enter the convolvers at t = 0, the 
unknown phase position will overlap the correct phase re- 
ference of one of the N convolvers at "some time during the 
interval f£ [MT Ct {2M — l)T c ] = [T, 2T — T c ] . After sam- 
pling, the largest of the resulting MN samples is chosen as the 
correct phase of the incoming waveform. Since one of the 
MN samples is the correct sample, it is therefore possible to 
initially acquire in 2T seconds. The actual mechanism for 
making the final decision regarding acquisition will be de- 
scribed in the next section. 

Suppose now that the initial phase uncertainty K is greater 
than MN. For this case, the recirculating delay line, shown in 
Fig. 1 , is used at the input to the system with both switches 
A and B closed. On. the input side, after the first T-second 
segment of the received waveform is in the system, switch A 
is opened and the 7"-second segment is recirculated / times, 
where J equals [K/MN] , and [x] is the smallest integer greater 
than or equal to x. On the reference side, the first N subcodes 
are put through the respective convolvers during the inter- 
val [0, TJ. From [7, 27], no reference input is used, but 
from 127, 37], the next N subcodes are used as reference 
inputs. In general, the /th set of N subcodes, * — I, J, is 
used as reference inputs during the interval [(/ - 2)7*, (/ — 1)7*] . 
For simplicity, in most of what follows, K is taken equal to 
MN, and hence 7=1. 

Considering next the lock mode, switch A is kept closed, 
switch B is kept open, and all other switches that were in 
position 1 are now put in position 2, so that each convolver 
output passes through the appropriate delay line. (Note that 
the recirculating delay line at the convolver input is not used.) 
In this mode of operation, the N outputs are noncoherent^ 
summed to yield an effective increase in integration time over 
that used in the search mode. To reset the reference input 
during the lock mode simply requires separating successive 
subcodes into any of the N convolvers by NT seconds. 

Notice that when the system is locked to the correct phase 
position, all the inputs to the delay lines have been despread, 
and hence, the delay lines themselves do not have to be broad- 
band devices. Similarly, if an incorrect phase position is being 
observed, while the inputs to the delay lines will no longer 
be despread, using narrow-band delay lines will still not re- 
sult in any performance degradation, since under this latter 
condition, a signal component out of the delay lines is not 
wanted. 
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IIL Performance Analysis 

Let the reference irTfjut to the zth convolver be denoted 
by PN ( ;_ 1)M+1 (r), where PN k (i) denotes a segment of a 
longer pseudonoise sequence [denoted PN (r)] from (Jt - l)r 
to (k - 1 + M)T Ct M is the-number of chips in each subcode*) 
• and T c is the duration of each chip. Therefore, PN (/ „ 1)M + i (r) 
is a subcode in the interval [(z - 1)MT C , iMT c \ t i =1,2, 
N. For simplicity in notation in what follows, the reference 
subcode PN (/ _ 1>M+1 (r) wiU be denoted by PN/*> (r). 
(Note that we have artibrarily defined the start of subcode 
PN x (*) (r) to correspond to t = 0.) 

To determine the probability that the correct phase posi- 
tion has been found, assume that PN,- (f) is the subcode 
that is currently being received, arid to be" specific, assume the 
received subcode PN/ (0 plus noise enters .the left side of the 
convolvers of Fig. 1 at .the same instant of time that the 
reference subcodes PN/*) (f - f). / = 1, •» Ar enter the right 
side. (Note that by making the convolver length equal to 2T, 
it is not necessary to have synchronism between the start of 
a received subcode and the start of the reference codes; it 
is done here only for convenience.) Also, for simplicity, as- 
sume J = 1. If we concentrate on the zth convolver, then from 
[1], it is straightforward to show that at any time r e [r, 
27*], the convolver output is given by 



if T x is set equal to 2T, then (1) reduces to 

: cos (2w 0 t - <p)+f PN,<*>(-2r + 3T + r) 
' °os (w 0 r)h(r) dr cos (2w 0 f - 0) 

+ / PN/ /? )(^2r+ 3!T+t) 

' sin 0v O 7-)n(r) dr sin (2w 0 r - 0). (2) 

Assuming initiaUy that z = / and that the sampling time T 
corresponds to precisely sampling the correct convolver at 
the peak of its response, then 



U;(T S )=[±AT+N C/ (T S )} cos [2 w 0 T s - 0j 
+ M S/ (T S ) $m [2w 0 T s - <t>] 



where 



PNy(r) cos w Q r + n(r)] 

PNfC«>(-2f + T l + r + r) 
cos [vv 0 (r - 2f) + 0] dr 



0) 



/T 
PN y ^)(r)n(T) cos w 0 r dr 

d 



i sin iv 0 r <77\ 



(3) 



(4) 



(5) 



where T A is the length of the convolver, X is a constant Both N C .(T S ) and N*XT.) are zero-mean r™««^ a 
amplitude, 0 is a random phase uniformly distributed in ^^Jtl ^ "jr jLZ J ™ m6an Gaussian rand <>™ 



amplitude, 0 is a random phase uniformly distributed in 
[0, 2tt], arid n(t) is AWGN of two-sided spectral density 
T? 0 /2. The specific limits of integration of (1) depend upon 
the degree of overlap between the received waveform and 
the reference code. If we assume the received waveform 
has completely occupied the convolver when the subcodes 
first enter, then for t <E [ r, 2T] , the iimits are 0 to 7*. 

If double frequency terms are ignored in the integrals and 



variables with variance r) 0 T/4. They are also statistically in- 
dependent. 

For i ¥= /, the output of the convolver at the sampling 
instant is given by 6 



Vi(T s ) = [±AR f/ + N C \(T S )) cos [2w 0 T s - <t>] 
+ ^.(rjsin f2 Wo ^-0] 



(6) 
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where N C .(T S ) and N S£ (T S ) have the same statistical properties 
as N CJ (T S ) and N S .(T S ), and where 



(7) 



MM- 1 

< 2 Hst>s Q }. 



If ^ £i>{S/ > 5y}, then P /y can be shown to be given by [2] 

=4 [1 - QCVb,,^) + QCVSTx/F)! (9) 

where 



2 L . (2a*)* 

± [4a"(l-p^)]'/2j 



(10) 



(11) 



(12) 
(13) 



and /oC*) is t« e modified Bessel function of the first kind 
and zeroth order. 

As a special case, if all the p ( y are approximately zero, 
then reduces to 



and (8) reduces to 



(14) 
(IS) 



In (14) and (1 5), E = A 2 T/2 and is the energy per integration 
interval, of the transmitted signal. As a final observation 
for the search mode analysis, note that for this special case 
(all = 0), the exact probability" of error can be found [3] 
and is given by 



Finally, with an appropriate change of notation, (6) can be 
seen to apply to the output of any of the N convolvers at 
any of the M — I incorrect sampling times. Therefore, if the 
output of the (th square-law detector at the kth sampling 
time is denoted by S^T + kT c ) p and if the correct sampling 
time (which equals 3T/2 for the specific timing being used 
here) corresponds to k — k c> then an error will occur if 

Si(T + kT c ) > S f (T + k c T c ) for any k, i * j, 
or if 

S / (T±kT c )>S / (T+ k c T c ) for Jfc * k c . 

These conditions correspond to an incorrectly matched con- 
volver at any sampling time having an output larger than that 
of the correct convolver at the correct sampling time, and to 
the correct convolver having a larger output at some time other 
than the correct sampling time, respectively. 

Denoting by {S m }, m - 1, *•*, MN - 1, the set of samples 
of either S&T + kT c ) for any k or S^T + kT c ) for k * k Ct 
and denoting Sj{T + k c T c ) by 5 0 , then an error occurs if 
$rn > So anv m ' Hence, the probability of error in the 
search mode, P e ^, can be upper bounded by 

= P{(Si > S 0 )or (S 2 > S 0 ) or »• or (S M „_ , > S 0 )} 



-'-TC"."') 5 

T E_ MN - i - \ ~\ 
L""no MM - i J 



M — 1 -i 



NM - 1 



(16) 



Consider now the lock mode, in which all the switches 
in Fig. 1 are set at position 2. In this latter mode, the N 
convolver outputs are all appropriately delayed and summed, 
but only at the specific phase position decided upon in the 
search mode. 

To analyze the performance of the system in the lock 
mode, the approximation that all the p,y = 0 is made in 
order to simplify the analysis. With this approximation, it 
can be seen that when an incorrect phase position is being 
observed, the decision variable has a central chi 2 distribution 
with 2N degrees of freedom, and when the correct phase 
position is being observed, the decision variable has a noncen- 
tral chi 2 distribution with 2N degrees of freedom and noncen- 
trality parameter a x 2 = NA^T*. Therefore, the probabilities 
of false alarm and correct detection in the lock mode are 
given by 



(2v 0 TYT(N) 



yM-l € -X/2-n 0 T dy 



(17) 



(8) and 



M) 



e ~(y+a 1 2 ) /2 T)or 



(18) 



respectively, where 7 is the threshold, T(-) is the gamma 
function, and /jv_ x (-) is a modified Bessel function. 

IV. Search Mode and Lock Mode Transitions 
In order to complete both the system description and the 
performance analysis, the specific procedure for transition- 
ing from one phase position to the next in the search mode 
must be specified, as well as the procedure for transitioning 
from the search mode to the lock mode. Towards this end, 
a search-lock strategy similar to that described in [4] will 
be used. In particular, once a phase position has been decided 
upon in the search mode, two more successive **hits** at that 
phase position will be required in order to enter into the lock 
mode. That is, once ail phase positions within the uncertainty 
range have been examined and an initial decision made as to 
which is the correct position, the procedure is repeated two 
more times. If the same relative phase position is chosen 
both times, the system advances to the lock mode. Similarly, 
two successive "misses 0 in the lock mode will be necessary 
in order for that phase position to be rejected and the search 
mode reinitiated. However, in the lock mode, only the specific 
phase position previously decided upon in the search mode is 
examined, not all possible phase positions falling within the 
uncertainty region. 

A flow diagram of this procedure is shown in Fig. 2. It is 
well known (see, e.g., [4], [5 J) that this type of process can 
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12 3 1 5 6 



Fig. 2. Markov chain model of search/lock strategy. 

be modeled as a finite. Markov chain with absorbing bound- 
aries. From Fig. 2, it is seen that there are six states to this 
particular Markov chain, labeled 1-6- States 1 and 6 are 
designated R for "rejection/* meaning when either state is 
reached, the phase position currently being examined is re- 
jected (i.e., the Markov chain has reached one of its absorb- 
ing states). States S x and S% correspond to the two search 
mode states, whereas LI and L2 correspond to the two lock 
mode states. 

In Fig. 2, P s denotes either P E W or 1 - P e M (depending 
upon whether the phase under consideration is incorrect or 
correct, respectively) and P^ represents Pf^ or P^W under 
similar circumstances for the lock mode. Using the absorbing 
Markov chain model for this system, the canonical form for 
the transition matrix, described in [4] and [5], is shown 
below. 
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1 -p* 


0 


Pi 


0 


4 


0 


0 


0 


.0 




1 -Pl 


5 


1 0 


1 -?L 


0 


0 


Pl 


0 



Forming the matrix 
1 



(19) 



I-Q 



(1-^) 
0 
0 



-Ps 
1 

0 
0 



-Ps 

-Pl 



0 
0 

-d -P L ) 
1 



(20) 



L = 1023 




I0" 6 

11.0 U.8 12.6 13.4 K.2 15.0 15.6 16.6 

Fig. 3. Upper bound on probability of error in the search mode for N - 11 
and M = 93. 

and [51, it is shown that if ry is the mean time to arrive at 
any absorbing state given an initial transient state /' (e.g., 
Si), then the vector of elements T/, denoted r, is given by 



T2 

^3 
U 



(22) 



where X is a vector whose elements are the times needed to 
make each transition. 

Since the time needed to move from one state to the next 
state for arbitrary J is 2TJ in the search mode and (N + l)T 
in the lock mode, it can be shown in a straightforward manner 
that the mean dwell time is given by 

T2= — (1-^ 11 -WW (23) 



and the expected time to loss of lock is given by 
2 -P L 

U = ~ — 0* + Or. (24) 



(1 -Pl) 2 



where / is an identity matrix, and defining the new matrix 

B— [I - Q]~ lW R, 

it is shown in [4] that the probability of entering lock is given 
by entry b 2 6 of tlie above B matrix. Notice that the notation 
b^26 refers to elements of the R matrix.- That is, the rows of 
B are numbered 2, 3, 4, and 5, and the columns of B are num- 
bered 1 and 6. For the system being considered 1 here, the 
probability of entering lock is given by 



*>2 6 — " 



Ps 2 



1 -Ps) 



(21) 



In addition to the probability of entering lock, other quan- 
tities- of interest are the mean dwell time in an incorrect phase 
position and the expected time to loss of lock, given that the 
system is in state L2 at the correct phase position. From [4] 



V. Numerical Results and Discussion 

As a specific example, it will be assumed that a PN se- 
quence of length L = 1 023 is being received simultaneously 
in N = 1 1 convolvers, thereby resulting in M = 93 unknown 
phase positions in each convolver. If the approximation 
leading to (14M16) is used, namely, that all the p if of (13) 
are zero, the curve of P e ^ s ) shown in Fig. 3 results. For the 
same system in the lock mode, Fig. 4 shows results of P^ 
for various values of 7 (or equivalently for various values of 

The variations of P f o and P d with iV* are illustrated in Figs. 
5-9. Fig. 5 shows P fa W versus Efa 0 for three values of N t 
N = 5, 10, and 15. Similar results for P d (0 are presented in 
Fig. 6, and Figs. 7-9 show how P d (0 varies with P f& (0 for a 
given value of N> 

As a perspective on the effect of approximating the p,y 
by zero for i ^ /, the exact P if of (9) were computed, and some 



BNSDOCID: <XP 22941 15A_L> 



MIJLSTEIN era/.: RAPID ACQUISITION FOR DIRECT SEQUENCE SPREAD-SPECTRUM 



597 




del 



3 4 5 6 

E/77 0 (dB) 

Fig. 4. Probability of detection versus B/rj 0 for various values of threshold, 
for N - 1 1 and M = 93 in the lock mode. 
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Fig. 7. Probability of detection for various thresholds in lock mode (N 
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Fig. 9. Probability of detection for various thresholds in lock mode (N = 
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Fig. 10. Probability of error for various relative starting phases of reference 
and input for N = 1 1 and M = 93 in the search mode. 



of the results are shown in Fig. 10. For the same PN sequence 
of length L — 1023 referred to above, a specific subsequence 
of length M = 93 was fed into all 1 1 convolvers and the re- 
sulting .10*23 values of p /; - were computed. That is, the starting 
phase of the reference input sequence PN/*) (f) took values 
0, 937* c , — , 930T C , corresponding to (i = 1, 2, U), where i 
is the index of the convolver. Also, the phases of PN; (r), 
the input sequence, took values 0, T c> 2T C , — , 92T C , cor- 
responding to / = 1, 2, 3, ***, 93. Hence, p (/ - was calculated 
for all these combinations of / and ; over the observation 
period T, and these p,- ; were then used to compute the cor- 
responding values of P^. 

There are four curves of along with their corresponding 
values for p l7 shown in Fig. 10. These curves correspond to 
0, /) pairs of (2, 1), (3, 1), (4, 1), and(i, 3). With the notation 
introduced above, (2, 1) corresponds to the P^ of the second 
convolver when the starting phase of its input signal is zero 
and that of its reference is 937V Similarly, (1, 3) corresponds 
to the first convolver when its input has a starting phase of 
2T C and its reference a starting phase of zero. The maximum 
variations in P if - t shown by these curves at any given SNR, 
come very close to the largest variation in P t j when all pos- 
sible (/,/) combinations are considered. 

Also plotted on this figure is the curve corresponding 
to Pij = 0 (i.e., I e -(i/2)SNR) Upon comparing these curves, 
it can be seen that the probability of error generated by 
approximating the p ;/ - by zero falls below all the other curves 
and, indeed, can be shown analytically to be a lower bound 
on performance. Similarly, the worst-case performance can be 
shown to correspond to that set of subcodes whose cross 
correlation is maximum in absolute value. For the spreading 
sequence that we chose, curve "tf" on Fig. 10 yields worst- 
case results. 

The advantage of the scheme, of course, is that in IT 
seconds, MN phase positions are examined. In a standard 
serial* search technique, examining MN phase positions re- 
quires MNT seconds, and so an improvement of a factor of 
MN/2 is achieved. To be more specific, consider a serial 



search scheme that accumulates all the correlation samples 
over the entire uncertainty range and then makes a decision 
on the basis of the largest sample. If the integration time for 
this latter procedure is T seconds per phase position in the 
search mode and NT seconds per phase position in the lock 
mode, then, assuming the search-lock strategy of Fig. 2 is 
employed, the mean dwell time is given by 



T2 (l-^) 2 (l T 



(25) 



where, as before, K is the total phase uncertainty. 

For either the parallel convolver scheme or the serial search 
scheme described above, the average time to correct acquisition 
can be found as follows. Assume that the time spent in lock 
at an incorrect phase position each time: a false lock occurs 
is given by T ff where T f is" given by (23) for the parallel convol- 
ver scheme and (25) for the serial search technique. Also, let f 
be the time spent making a decision at any step in the search 
mode. Note that T - KT for the serial search and f = 2TioT 
the parallel scheme, assuming J = 1. Then the mean time to 
acquire is given by 

T acq = f(l - P e &) + (T f + 2T )(1 - P e W)P e i*) 
+ (2 T f + 3f)(l - /> e (*>)(/> e C*>)2 
+ -+ [«?>+(« + l)f](l -P e &) 
* [/' e ( j) ] m + *" 



_ T+TfP e (*) 
1 - P e M 



(26) 



BNSDOCID: <XP 



_2294115A_I_> 



V. 



4 

1 



MILSTEIN etal: RAPID ACQUISITION FOR DIRECT SEQUENCE SPREAD-SPECTRUM 




Fig. 1 1 . Mean dwell lime versus the false alarm probability in the search 
mode for various values of the false alarm probability in the lock mode. 
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Fig. 12. Mean time to loss of lock versus the probability of correct detection 
in the lock mode. 
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Notice that if T f P e W < f t then 

2 

r acq (parallel convolver) ^ — r acq (serial search). (27) 

It is, of course, necessary to have an acceptable probability 
of error in the search mode, and if the E/r} 0 required to yield 
a given level of performance is too large, one can always in- 
crease the number of states in the search mode. While this will 
not decrease P e ^ on any given try, it will decrease the prob- 
ability of entering lock at an incorrect phase position. Also, 
this further illustrates a point made in [4], namely, that a 
large integration time in the lock mode is desirable (recall 
that in the scheme described here, the integration time in the 
lock mode of any one convolver is the same as it is in the 
search mode, namely T seconds, but that noncoherently 
summing N output samples has the effect of yielding a larger 
integration time). Since one has to search all the phase posi- 
tions in the uncertainty region in the search mode, but one 
only has to check a single phase position in the lock mode, 
once one has entered the lock mode, it is desirable to stay 
in that mode (i.e., not lose lock) for as long as possible. As 
long as the overall mechanism of entering the lock mode 
results in a correct decision with a reasonably high probabil- 
ity, the additional time spent in false lock caused by employ- 
ing this strategy can be kept to a tolerable level. 

Finally, it is at times more meaningful to consider such 
parameters as mean dwell time at an incorrect phase position 
and mean time to loss of lock at the correct phase position 
than it is to consider the probabilities of detection and false 
alarm. (Which set of parameters is more relevant to a given 
system is typically a function of how long the acquisition 
sequence is available, with the former parameters being more 
.significant in continuously running sequences, as in certain 
navigation systems, and the., latter parameters being more 
significant in more bursty systems, such as packet radios.) 
To illustrate the behavior of this particular acquisition scheme 
with respect to mean dwell time and mean time to loss of lock, 
consider Figs. 11 and 12. In Fig. 11, Tj/r is plotted versus 
the false alarm probability in the search mode, with the false 
alarm probability in the lock mode as a parameter. In Fig. 12, 
r 4 /r is plotted versus the probability of correct detection in 
the lock mode. 

VI. Conclusion 

A technique for the rapid acquisition of a OS spread- 
spectrum signal has been presented. The scheme relies on 
parallel processing different subcodes of a longer spreading 
sequence simultaneously in SAW convolvers, and results in 
a reduction of search time which is proportional to M times 
//, where M is the number of chips being processed in each 
convolver and N is the number of convolvers. Results were 
presented for such key system parameters as probability of 
detection and mean time to loss of lock. While only analytical 
results were considered in the paper, it is believed that the 
system could be implemented with current technology in 
SAW devices. 
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